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ABSTRACT 

Population III stars can regulate star formation in the primordial Universe in 
several ways. They can ionize nearby halos, and even if their ionizing photons are 
trapped by their own halos, their Lyman-Werner (LW) photons can still escape and 
destroy H 2 in other halos, preventing them from cooling and forming stars. LW escape 
fractions are thus a key parameter in cosmological simulations of early reionization 
and star formation but have not yet been parametrized for realistic halos by halo or 
stellar mass. To do so, we perform radiation hydrodynamical simulations of LW UV 
escape from 9-120 M 0 Pop III stars in 10 5 to 10 7 M 0 halos with ZEUS-MP. We 
find that photons in the LW lines (i.e. those responsible for destroying H 2 in nearby 
systems) have escape fractions ranging from 0% to 85%. No LW photons escape the 
most massive halo in our sample, even from the most massive star. Escape fractions 
for photons elsewhere in the 11.18-13.6 eV energy range, which can be redshifted into 
the LW lines at cosmological distances, are generally much higher, being above 60% 
for all but the least massive stars in the most massive halos. We find that shielding of 
H 2 by neutral hydrogen, which has been neglected in most studies to date, produces 
escape fractions that are up to a factor of three smaller than those predicted by H 2 
self-shielding alone. 

Key words: early universe - cosmic background radiation dark ages, reionization, 
first stars - stars: Population III - radiation: dynamics. 


1 INTRODUCTION 

The first stars form at 2 ~ 20-30, or about 200 Myr after the 
big bang. Early numerical simulations suggested that Pop III 
stars are very massive, 100 - 500 M 0 , and form in isolation, 


one per halo (Bromm, Coppi & Larson 

1999 2002 

Abel, 

Bryan & Norman(2000 2002, Nakamura & Umemura 

2001 

Glover 2005 

. More recent models indicate that most Pop III 


stars form in binaries (Turk, Abel & O’Shea 2009) 0 

in small 

clusters ( 

Stacy, Greif & Bromm 2010 Clark et al. 

2011a,b; 

Greif et al. 2011 2012 Smith et al. 2011 Dopcke et al. 2013; 

Glover 

2013 

Stacy & Bromm 

2014 

1 with a wide range of 


masses, all the way from the sub-solar regime to ~10^ M 0 . 
These calculations suggest that the final mass of a star is 


either limited by dynamical processes (Clark et al. 2011a; 


Greif et al.[2012 

Smith et al.|2012) or by radiative feedback 

(McKee & Tan 

|2008) Hosokawa et al. 

2011 2012 Stacy, 

Greif & Bromm 

2012; Hirano et al.12014). 


If any Pop III stars formed with masses below ~0.8 M 0 , 


E-maibschauer@uni-heidelberg.de 


they will have survived until the present day, and they may 
therefore be detectable in surveys targeting extremely metal- 
poor stars in the Galactic bulge and halo as well as in nearby 


satellite galaxies (see the reviews by Beers & Christlieb 2005 
Frebel|2010 1. This could potentially allow us to set limits 


on the low mass end of the primordial stellar initial mass 
function (IMF, see e.g. Tumlinson 2006 Salvadori, Schnei- 


der fe Ferrara|20~07 jHartwig et al.|2015 1. Unfortunately, the 

direct detection of Pop III stars in the high redshift Universe 
is unlikely even with future 30 meter telescopes (although 
Rydberg, Zackrisson & Scott 2010). As massive stars 


are very short lived, the upper end of the Pop III IMF is 
thus more difficult to constrain. This can only be done in¬ 
directly, for example by observing and analyzing the super¬ 
novae (SNe) that mark the end of the lives of very massive 
stars. These may indeed be visible to the James Webb Space 
Telescope ( JWST ), the Thirty-Meter Telescope (TMT), the 
Extremely Large Telescope (ELT) and the Giant Magellan 
Telescope (GMT) (e.g., Whalen et al.|[2013a|b| 2014| Smidt 
et al. J2°14|. Another promising possibility is again Galactic 
archeology. The comparison of the nucleosynthetic yields of 
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fe Bromm|2008 Joggerst et al.|2010l. 


Pop III SNe to the chemical abundances of extremely metal- 

the star (Susa 

2007 

Whalen et al. | 2008 Whalen, Hueck- 

poor stars observed in our Milky Way suggests that many 

Pop III stars may had masses in the range 15 — 40 Mq (e.g., 

staedt & McC 

onkic 

2010;; see also Shapiro, Iliev & Raga 

2004} Iliev, Shapiro & Raga||2005. Hasegawa, Umemura & 

Beers & Christliebi,2005 Frebel et al. 2005} |Frebel, Johnson 

Susa|2009; Susa, Umemura & Hasegawa, 2009). 


In this paper, we focus on the high-mass end of the 
Pop III IMF. Massive Pop III stars had the capability to 
profoundly transform their environment due to the strong 
radiative, chemical, and mechanical feedback they provide. 
They can evaporate the gas from their dark matter (DM) ha¬ 
los, engulf nearby halos with both ionizing and LW band UV 


radiation l 

Whalen, Abel & Norman|2004 

Alvarez, Bromm 

& Shapiro 

]2006; Abel, Wise & Bryan||2007j Whalen et al. 

2008 

), and alter their composition with the metal they pro- 


duced and expelled. LW photons lie in a series of discrete 
lines located in the energy range 11.18 - 13.6 eV (corre¬ 
sponding to a wavelength range of 1110 A - 912 A ) and can 
travel great distances in the primeval universe because they 
lie below the ionization threshold of atomic hydrogen. They 
can delay or suppress star formation in nearby halos because 
they photodissociate H 2 and prevent the gas from cooling 


(see e.g. Haiman, Rees & Loeb 1997; Haiman, Abel & Rees 

2000 

Glover & Brand 2001, 2003; |Machacek, Bryan & Abel 

2001 

Susa & Umemura| 2006, Susa 2007| O’Shea & Nor- 

man 

2007 Wise & Abel 

2008). How Pop III stars regulate 


subsequent star formation is central to the rise of stellar pop- 


ulations in the first galaxies iCiardi & Ferrara||2005 

John- 

son et al.||2009, Greif et al.||2010 Jeon et al.|2012| 

^wlik, 

Milosavljevic & Bromm|2013 

Wise et al. 20121. This begins 


with understanding how the strength of the LW background 
changes over time, which depends on the masses of Pop III 
stars and the halos in which they reside. 

Over the past two decades, many studies have examined 
the influence of LW and ionizing radiation on the formation 
of stars at high redshift. They found that a photodissoci- 
ating background cannot halt Pop III star formation but 
can delay it (e.g., Machacek, Bryan fe Abel] 2001 |Wise fe 


Abel|2007| |Q’Shea fe Norman|2007[ |Safranek-Shrader et al. 

2012b. An early X-ray background created by accretion onto 


stellar remnants and by emission from supernovae can cat¬ 
alyze additional H 2 formation that offsets the effects of LW 
radiation to some extent, but radiative feedback overall is 


dominated by the LW background (Glover & Brand 2003 
Machacek, Bryan & Abel 20031. The evolution of this back¬ 


ground over time in general seems to play only a mi nor r ole 
if it changes less rapidly with redshift a than 10 -2 ^ 5 (Visbal 


et al.| 2014). Simulations suggest that Pop III stars forming 
in halos illuminated by significant LW backgrounds will be 
more massive because their host halos must grow to larger 
masses before beginning to cool (O’Shea & Norman 2007 


Hirano et al.|[20l5 Latif et al.||2014 |. When they do begin 

to cool, collapse rates at their centers are higher and lead to 
more massive stars. 

The effect of a single Pop III star on a nearby halo 
has been examined in much greater detail than in large cos¬ 
mological boxes. Radiation hydro dynamical simulations by 
Susa & Umemura ( 2006| ) find that star formation in the 
vicinity of a source star is possible if the star forming region 
exceeds a density threshold (see also Glover fe Brand|2001 1. 
Other studies in this vein show that star formation can be 
either promoted or suppressed depending on the mass of the 
halo, the mass of the star, and the proximity of the halo to 


In some cases the H 11 region of a Pop III star may fail to 
break out from a halo because of large central gas densities. 
In such cases, LW photons might still exit the halo because 
their energies lie below the ionization limit of hydrogen, even 
if H 2 self-shields against this flux to some degree. Most sim¬ 
ulations of radiative feedback with Pop III stars in cosmo¬ 
logical boxes simply assume uniform LW backgrounds. A 
parametrization of LW escape fractions from primordial ha¬ 
los as a function of halo and stellar mass could provide the 
strength of this background from first principles in future 
simulations. However, almost no work has been done to pro¬ 
duce such a parametrization. In the only previous study on 


this topic of which we are aware, Kitayama et al. (2004) ex¬ 


amined LW escape fractions from Pop III star-forming halos 


with a Navarro, Frenk & White (19971 radial density profile, 


varying the mass of the halo and central star. 

In an effort to improve on this previous study, we per¬ 
form radiation hydrodynamical calculations of LW escape 
from Pop III minihalos with the ZEUS-MP code. We then 
post-process the simulations with semi-analytical methods 
and calculate escape fractions in two limits, the near-field 
and the far-held. In the near-Held case - the only scenario 


considered by Kitayama et al. (2004 1 - we compute the es¬ 
cape fraction of photons in the LW lines themselves. This is 
the value that is relevant if we are interested in the effect 
of radiation from the halo on H 2 in its immediate vicinity. 
In the far-held case, on the other hand, we are interested in 
the total fraction of the photons lying between 13.6 eV and 
11.2 eV that can escape from the halo. This is the important 
quantity if we are interested in the effect of radiation from 
the halo on H 2 located at cosmological distances. 

In both cases, we account not only for H 2 self-shielding 
but also for shielding from the Lyman series lines of atomic 
hydrogen, which were not considered by |Kitayama et al \ 
12004). In some circumstances, this can significantly affect 
the LW escape fraction. 

Our paper is structured as follows. In Section 2 we de¬ 
scribe our numerical models. In Section 3 we tabulate LW 
escape fractions by halo and stellar masses and we conclude 
in Section 4. 


2 METHOD 

To calculate LW escape fractions for a given star and halo 
we hrst evolve the H 11 region of the star with the radiation 
hydrodynamics simulation code ZEUS-MP. We then post¬ 
process the prohles for the ionization front (I-front) and the 
surrounding halo with semi-analytic calculations to deter¬ 
mine how many LW photons exit the halo. LW escape frac¬ 
tions in both the near-held and far-held approximations are 
considered, as described below. 


2.1 ZEUS-MP 

ZEUS-MP is an astrophysical radiation-hydrodynamics code 
that self-consistently couples photon-conserving raytracing 
UV transport and nonequilibrium primordial gas chemistry 
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to gas dynamics to evolve cosmological I-fronts (Whalen & 


Norman 2006 


2008a|b ). We evolve mass fractions for H, H + , 


He, He + , He 2+ , H“, H(T, H2, and e“ with nine additional 
continuity equations and the nonequilibrium rate equations 
of Anninos et al. (1997) in which the species are assumed to 
share a common velocity distribution. Mass and charge con¬ 
servation, which are not formally guaranteed by either the 
network or advection steps, are enforced at every update 
of the reaction network. Heating and cooling due to pho¬ 
toionization and chemistry are coupled to the gas energy 
density with an isochoric update that is operator-split from 
the fluid equations. Cooling due to collisional ionization and 
excitation of H and He, recombinations of H and He, inverse 
Compton scattering (IC) from the CMB, bremsstrahlung 
emission, and H2 cooling are all included in our models. The 
chemical and cooling rate coefficients are the same as those 
used in Anninos et al. (19971, with one important exception: 
we use case B rates to describe the recombination of hydro¬ 
gen and helium, rather than case A as used in Anninos et al.| 
( 1997 1. 

We use 120 energy bins in our photon-conserving UV 
transport scheme, 40 bins that are uniform in energy from 
0.755 to 13.6 eV and 80 bins that are logarithmically spaced 
from 13.6 eV to 90 eV. We normalize photon rates in each 
bin by the time-averaged ionizing photon rates and sur¬ 
face temperatures for Pop III stars from Tables 3 and 4 of 


Schaerer (20021. The radiative reactions in our simulations 


are listed in Table 1 of Whalen & Norman (2008a I, and the 
momentum imparted to the gas by ionizations of H and He 
is included in the photon transport. Photon conservation is 
not used to calculate H2 photodissociation rates. They are 
derived along radial rays from the star with the self-shielding 
functions of Draine & Bertoldi (1996 hereafter DB96) mod¬ 
ified for thermal broadening as a proxy for the effects of gas 
motion (equations 9 and 10 in Whalen & Norman 2008aI, 
with r~ 2 attenuation taken into account. 

Two-photon emission from recombining He + and He 2+ 
can produce photons capable of photo-dissociating H2. How¬ 
ever, we do not account for the effects of this nebular emis¬ 
sion, as it is unimportant in comparison to the effects of the 
direct stellar emission. 


2.2 Halo Models 


We adopt the one-dimensional (ID) spherically-averaged 
halo profiles used in |Whalen, Hueckstaedt fe McConkie| 
12010), which were based on the results of cosmological sim¬ 


ulations carri ed out wi th t he En zo adaptive mesh refinement 
(AMR) code (Bryan et al.|2014 1. Their masses are 6.9 x 10 5 
M 0 , 2.1 x 10° M 0 and 1.2 x 10 7 M 0 . This corresponds to 
the range of halo masses in which Pop III stars are expected 
to form via H2 cooling, and all three halos form at z ~ 20. 
Densities, velocities, temperatures and mass fractions for all 
nine primordial species are mapped from these profiles onto 
a ID spherical grid in ZEUS-MP. We summarize the proper¬ 
ties of these halos in Table |T] and plot their initial densities, 
velocities, temperatures, and species mass fractions in Fig¬ 
ure [T] 

The gas densities have nearly a power-law profile with 
slopes around -2.1, with the most massive halo having the 
highest densities. Temperatures vary from a few hundred 
K to several thousand K and generally increase with halo 


mass. This is to be expected since the virial temperature 
of a halo scales with the halo mass as TVir oc M 2 ^ 3 . Virial 
shocks are visible in all three profiles at ~ 100 pc, where 
infall velocities abruptly decrease as accretion flows crash 
into increasingly dense regions of the halo. Shock heating at 
these radii is also evident in the temperature profiles, and is 
almost strong enough to collisionally ionize atomic hydrogen 
in halo 3. 

ZEUS-MP does not evolve DM particles as in cosmolog¬ 
ical codes such as Enzo or GADGET; instead, an additional 
gravitational potential is implemented as a proxy for the 
DM potential of the halo. We take its potential to be that 
required to keep the baryons in hydrostatic equilibrium on 
the grid. The mass associated with this potential is nearly 
the same as that of the halo. We interpolate this precom¬ 
puted potential onto the grid at the beginning of the run 
but it does not evolve thereafter. Taking the potential to 
be static is a reasonable approximation because merger and 
accretion timescales at z ~ 20 are on the order of 20 Myr, 
significantly longer than the lifetimes of most of the stars 
that we consider in this study (see Table [2J. We therefore 
do not expect the dark matter distribution to evolve much 
over the lifetime of the star. 


2.3 Hu Region Simulation Setup 

We center the Pop III star at the origin of a ID spherically- 
symmetric grid that has 500 ratioed zones in radius. The 
inner boundaries are at 2.0 x 10 1 ' cm (halos 1 and 2) or 
4.0 x 10 17 cm (halo 3) and the outer boundaries are at 
slightly more than twice the virial radius of the halo. Re¬ 
flecting and outflow conditions are imposed on the inner 
and outer boundaries, respectively. Densities, energies, ve¬ 
locities and species mass fractions from our Enzo profiles are 
mapped onto the ZEUS-MP grid with a simple linear inter¬ 
polation of the logarithm of the given variable. We consider 
Pop III stars with masses 9, 15, 25, 40, 60, 80 and 120 M 0 
in all three halos, leading to a total of 21 models. The prop¬ 
erties of these stars are listed in Table 2 (|Sch aerer||2002| ). 
Each simulation is run out to the end of the life of the star. 

In two of the 21 simulations (the 15 M 0 star in halos 
1 and 2) the I-front is confined to small radii that require a 
finer grid in order to be resolved. We therefore use a grid of 
1000 zones with two contiguous blocks: a very finely spaced 
uniform grid with 600 cells to resolve the I-front followed 
by a ratioed grid with 400 cells that again extends to twice 
the virial radius of the halo. We ensure that the length of 
the innermost zone of the outer block is within 20% of that 
of the outermost zone of the inner block to avoid spurious 
reflections of shocks at the interface of the two blocks. 


2.4 LW Escape Fractions 

Molecular hydrogen has no dipole moment and is there¬ 
fore destroyed in a two-step photodissociation process, the 
Solomon-process (Stecher & Williams 19671: H2+7 —> H2 —» 
2H. Incident radiation can excite the molecule from the 
ground state to an electronically excited state. Two of these 
excited states, B 1 Eu (known as the Lyman state) and C 1 ^ 
(the Werner state), are separated from the electronic ground 
state by less than 13.6 eV, and transitions to these states 
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Figure 1. Initial densities (upper left), temperatures (upper right), radial velocities (lower left) and species mass fractions (lower right) 
for the three halos. 


can therefore be brought about by photons with energies 
below the Lyman limit of atomic hydrogen. Within the Ly¬ 
man and Werner states, a variety of bound rotational or 
vibrational levels are accessible, and so transitions from the 
ground state to the Lyman or Werner states occur through 
a series of discrete lines, known as the Lyman-Werner band 
system, or simply as the Lyman-Werner lines. Once excited 
by a photon in one of these lines (a LW photon), the H 2 
molecule remains in the Lyman or Werner state for only a 
very short time At ~ 1CP S s, before decaying back to the 
electronic ground state. Most of the time, the H 2 molecule 
decays to a bound ro-vibrational level in the ground state. 
However, a small fraction of the time, the decay occurs to 
the vibrational continuum, resulting in the dissociation of 
the H 2 molecule. The dissociation probability depends on 
the details of the incident spectrum and the density and 
temperature of the gas (which fix the initial rotational and 
vibrational level populations), but is typically around 15% 
(DB96). The remaining 85% of LW photon absorptions re¬ 
sult in decay back to a bound state. In a small fraction of 
cases, this results in a photon with the same energy as the 
original LW photon (albeit with a random direction). How¬ 
ever, the majority of the time, the re-emitted photon has too 


|Brand|[2001[ for more details, see the discussion in Section 
3.4 of). 

When the H 2 column density is large, the gas becomes 
optically thick in the LW lines. This can prevent radiation 
from an external source from reaching the center of a halo or 
LW flux from the center of a halo from escaping it. The gas 
therefore self-shields against LW radiation. In addition, LW 
photons can also be absorbed by most of the Lyman series 
lines of atomic hydrogen [^However, because the frequencies 
of the LW lines of H 2 do not coincide particularly closely 
with those of the Lyman series lines of H, this effect only be¬ 
comes important when the atomic hydrogen column density 
is very high, Nh ~ 10 2i cm -2 or above, so that the Lyman 


series lines are strongly Lorentz-broadened (Wolcott-Green 
fe Haiman||2011 1 . 


small an energy to bring about photodissociation (Glover & 


In the immediate vicinity of a halo, the photons re¬ 
sponsible for photodissociating H 2 are those emitted in the 
LW lines in the rest-frame of the halo. The absorption of 
these photons by H 2 and H within the halo can be conve¬ 
niently parameterized by a simple self-shielding function (see 
e.g. DB96; [Wolcott-Green fc Haiman|2011[ hereafter simply 


1 The exception is Lyman-a, which is located at too low an en¬ 
ergy. 
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halo i?200 (p c ) mass (Mg) 

1 256.8 6.9 x 10 5 

2 339.6 2.1 x 10 6 

3 495.2 1.2 x 10 7 

Table 1. Virial radii and total (dark matter + baryon) masses 
of the three halos. 


M (M 0 ) 

lifetime (Myr) 

log L (L 0 ) 

log T eff (K) 

9 

20.2 

3.709 

4.622 

15 

10.4 

4.324 

4.759 

25 

6.46 

4.890 

4.850 

40 

3.86 

5.420 

4.900 

60 

3.46 

5.715 

4.943 

80 

3.01 

5.947 

4.970 

120 

2.52 

6.243 

4.981 


WH11). The escape fraction of photons in this limit, which 
we term the near-field limit, is then simply given by the value 
of this self-shielding function evaluated at the virial radius 
of the halo. At larger, cosmological distances, the LW flux 
is redshifted by the expanding Universe or Doppler shifted 
due to relative velocities between halos. As a consequence, 
the LW lines in the source frame may not coincide with LW 
absorption lines in the local frame of a halo. A simple esti¬ 
mate of the escape fraction in this limit, which we term the 
far-held limit, is given by the fraction of the LW range 11.2 
- 13.6 eV not lying inside the combined equivalent width 
of the absorption lines. In the sections below, we describe 
in more detail how we compute the LW escape fractions in 
both limits. 


2-4-1 Near-Field Limit 


A LW photon escapes a halo if it does not photodissociate 
an H 2 molecule or is not absorbed and re-emitted as lower- 
energy photons. The LW escape fraction, f esc , can therefore 
be equated to the fraction of H 2 molecules that are shielded 
from either process. The factors by which H 2 is shielded 
from LW photons by other H 2 molecules and by H atoms 
are / a ^f eld and ie i d , respectively. Both factors account for 
all processes by which LW photons are absorbed, not just 
those that result in a photodissociation. The total factor by 
which H 2 is shielded can be taken to be the product of these 
two factors (see the commentary on equation 13 in WH11), 


/*H 2 rH 
J shield J shield • 


(i) 


To construct / s ^f eld we first determine the column density of 

H 2 , 


ri<i ?200 

Nh 2 = niXn 2 ,i(fi-ri-i), (2) 

2 = 1 

where rii is the number density of all particles and ri is 
the outer radius of the i th cell, where ro is zero. \H 2 ,i = 
A r H 2 ,i/A 7 totai,i is the H 2 particle fraction with Nu 2 ,i the num¬ 
ber of H 2 molecules and Wotai,* the total number of particles 
in bin i. The sum extends out to the virial radius of the halo, 


R200 = R{P = 200 f2b,0 Pcrit), (3) 


where p is the gas density and p cr i t = 3 H 2 (z = 20) / (8nG) ~ 
2.349 x 10 _26 g cnr J is the critical overdensity of the Uni¬ 
verse at z = 20. The virial radii are listed for all three halos 
in Table [T| 

The Doppler broadening of lines must be taken 
into account in /S eld . The thermal component of the 
Doppler broadening parameter, 6 d,t, is associated with the 
Maxwellian velocity distribution of the atoms and molecules, 


Pv (v)du 


I m 

v^^ exp 


/ mv 2 \ 

\ 2k bTJ ' 


(4) 


Table 2. Pop III stellar properties ||Schaerer|2002). 


From DB96, the Doppler broadening parameter is defined by 
b = FWHM/(41n2) 1/2 , where FWHM is the full-width half 
maximum of P v ( v) and is related to the standard deviation, 
a, of P v (v) by FWHM = 2y/2\E2a. This yields 


&D,T 


2 kgT 

mg 2 


(5) 


where kg is the Boltzmann constant and jtih 2 is the mass of 
an H 2 molecule. The temperature T in equation © is the 
H 2 -weighted mean temperature T e ff of the gas 


Teff = 


E ? = T 200 rn XH 3 ,iTiVi 

j:ZT 200 

EH!* 300 niXK 2 ,iTi%(r! -rU) 
EEl fl200 rn XH 2 ,i ^(rf - rfLr) ’ 


( 6 ) 


where T; is the temperature and Vi is the volume in cell i. 

The bulk motion of the gas can also affect whether or 
not a fluid element is shifted into the line of an outgoing LW 
photon. We account for this in an approximate way with an 
additional component to the Doppler broadening parameter, 


&D,T,v — Jb D T + bg — J —L -f 2cr 2 lt , (7) 

V y ?71 h 2 

where cr m t is taken to be the microturbulent velocity dis¬ 
persion, which is approximated as the H 2 -weighted velocity 
dispersion of the halo in analogy to T e a. We then calculate 
/shield w ith the fitting function in WH11 

,h 2 0.9379 

/shield ( H 2 ,bo,T, v ) (l + x / 1.879 

+ ( i° + °^o. 47 3 x exp[—2.293 x 10~ 4 Vl + x}, (8) 

where x = Ah 2 /(8.465 x 10 13 cm -2 ) and Dg 2 = 

&D,T,v/(10 5 cm s _1 ). 

We construct i e i d from the fitting function in WH11, 

/shield(A^ h) = , u 62 exp(-0.149 Xg ), (9) 

(1 -t- 

where ih = IVh/(2.85 x 10 23 cm -2 ). Escape fractions in the 
near-field limit are calculated for each profile of the Hu re¬ 
gion and surrounding halo as they evolve throughout the 
run, typically 1000 times during the simulation. 


2-4-2 Far-Field Limit 

In the far-held approximation, we begin by calculating the 
dimensionless equivalent width of all the LW lines, W'. It 
cannot exceed the total dimensionless width of the LW range 
11.2 - 13.6 eV: !F r ' rulx = ln(1110/912) ~ 0.2. The escape 
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Figure 2. Evolution of the I-front of the 60 Mg star in halo 1. The times are 0.175 Myr (black solid line), 0.35 Myr (blue dotted line), 
0.875 Myr (green dashed line) and 3.5 Myr, the end of the life of the star (red dash-dotted line). 


fraction in the far-field limit is then one minus the fraction 
of these two quantities, / esc = 1 — W'/W^x.- 

The transition from a lower level l (with vibrational 
quantum number v and rotational quantum number J) to 
an upper level u (with v' and J ', respectively) depends on 
the photoabsorption cross-section a u i and column density 
Ni. In a first step, we show the calculation for a single line. 
Later, we combine them to obtain the total equivalent width. 
We follow the notation of DB96. The equivalent width of the 
l to u transition can be written as 

WuiiNt) = Al - eM~Nia u i)]du. (10) 

If Nia„i < 1 we can treat the line width in the weak-line 
limit, in which exp(— Nia u i) ~ 1 — Nicr u i. The photodissoci¬ 
ation cross section can also be written in terms of the cross 
section at line center, a u i(v o), and a frequency dependent 
line profile function that is normalized to one. Therefore, we 
have 

W u i,w(Ni) ~ Ni(j u i{vq). (11) 

The shape of the equivalent width needs to be accounted 
for in the general case. Two processes determine the shape, 
Lorentz and Doppler broadening. 

Lorentz broadening, also known as natural broadening, 


is intrinsic to the line and is a consequence of the Uncer¬ 
tainty Principle: a quantum state with a lifetime At has 
an uncertainty associated with its energy A E such that 
AEAt ^ h/2. Consequently, radiative transitions to or from 
this state do not have a precise energy but instead occur with 
a range of energies with a distribution with a width ~ A E. 
Doppler broadening, also known as thermal broadening, is a 
consequence of the thermal motion of molecules and atoms. 
The Doppler shifts associated with this motion again lead to 
a spread in the frequency of the transition in the laboratory 
frame. If the gas particles have a Maxwell-Boltzmann veloc¬ 
ity distribution, the line profile is a simple Gaussian when 
Doppler broadening dominates. 

Doppler broadening usually dominates near the center 
of the line uq. However, the Doppler line profile function falls 
off exponentially away from the line center while the Lorentz 
profile falls off only as (v— vo)~ 2 . Lorentz broadening there¬ 
fore always dominates far enough from the centre of the line. 
Lorentz broadening can often be ignored far from the centre 
of the line because the line is so weak there that the optical 
depth due to that region of the line profile is -C 1. However, 
if Ni is very large this may no longer be true and we need 
to account for both Lorentz and Doppler broadening. 

When both effects are important, the line profile be- 
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comes a convolution of the Lorentz and Doppler profiles 
known as the Voigt profile. Instead of calculating the 
Voigt profile directly, we follow the approach of |Rodgers] 


& Williams (1974 hereafter RW74) and approximate the 


equivalent width of a single line to be 


W ul = [Wi >L + w, 


ul ,D - {W uU L * Wul,v!Wul,W?] 1 ' 2 (12) 


(equation 3 of RW74). Here, Wl is the equivalent width 
of the line, assuming only Lorentz broadening, Wd is the 
equivalent width for only the Doppler effect taken into ac¬ 
count, while Ww is the equivalent width in the weak line 
limit (equation 10 ). 

We use the approximation in RW74 to calculate Wd, 


Wd = A vdD(z), (13) 

where Ai/d = &d,t,vZ'o/c includes the Doppler broadening 
parameter 6d,t,v, as defined in Equation ( 6 ). The function 
D is a seventh-order polynomial given in the Appendix of 
RW74 and z = a u i{i/ 0 )Ni/{AuDn 1/2 ). 

We take the expression for Wl given in Belafhal ([2000), 
which is more accurate than the one from RW74, 


W L = 2 ttA ulL{z), (14) 


where Aul = T/2 is the Lorentz half-width of the line with 
the total de-excitation rate r, z = a u i{vo) Ni / {2 hAvl), and 
L is the Ladenburg-Reiche function used by Belafhal (2000). 
To combine the individual lines we first find their dimen¬ 
sionless equivalent widths. As the width of the line is small 
compared to its frequency, it can be written as 

W ul 

vo 


WL = 


(15) 


We assume an ortho-to-para ratio of 3:1 for molecular hy¬ 
drogen, so 7V or tho = 0.75 x Vh 2 and V para = 0.25 x Nh 2 . Ap¬ 
plying the transition rules, only one upper rotational level 
can be reached: J f — 1. We further limit our calculations 
to transition energies below the ionization limit of hydro¬ 
gen, E ^ 13.6 eV, because we assume that any photon with 
E ^ 13.6 eV ionizes a hydrogen atom. 

For each transition, we take the molecular data required 
to compute the dimensionless equivalent width - the oscil¬ 
lator strength and frequency of the transition, and the total 
radiative de-excitation rate of the excited state - from the 
papers by |Abgrall fe Roueff| ( |1989| ) and |Abgrall et alT] ( |1 992 [ ) . 

In the case of no line overlap, the total dimensionless 
equivalent width can be calculated by a sum over all indi¬ 
vidual lines, 


w ’ = J2 T, w L- ( 16 ) 

l u 


If there is overlap we account for it in the same way as DB96. 
To ensure that the total dimensionless equivalent width of 
our set of lines satisfies this constraint, we write it as 

W' = WL* [1 - exp(-W7W' lax )] . (17) 

Therefore, since f esc = 1 - / a bs = 1 - W'/W„ lax , we have 
/ esc = exp {-W'/WL^) • (18) 

As in the near-field case, we consider shielding by neu¬ 
tral hydrogen in addition to self shielding by H 2 . All Lyman 
lines starting from Ly,0 fall into the LW range and can there¬ 
fore reduce the LW escape fraction. Like WH11, we consider 


Lyman transitions (l = 1) from u = 2 up to u = 10. We cal- 

(20091 in analogy 


culate Wtr with data from Wise & Fuhr 


to W'. The two dimensionless equivalent widths are then 
summed and the total escape fraction becomes 

W' + Wh' 


/esc,H 2 &H = exp — 


W' 

v y m 


(19) 


3 RESULTS 

We first examine the evolution of the I-front and H 11 region 
in a halo and then present our results for / esc in the near¬ 
field and far-held limits. 


3.1 Evolution of the Hu Region 

The evolution of the I-front and H 11 region depends on the 
mass of the star and its host halo. At early times the I-front 
propagates very rapidly, leaving the gas behind it essentially 
undisturbed (an R-type front). It then decelerates as it ap¬ 
proaches the Stromgren radius of the halo. If the halo is 
massive and the star is not very luminous the I-front stalls 
at the Stromgren radius and advances no further. If the star 
is bright, the front may briefly loiter at the Stromgren ra¬ 
dius, but it then resumes its expansion driving a shock in 
front of it (a D-type front). As it descends the steep density 
gradient of the halo, the front can break through the shell 
and revert to R-type, flash ionizing the halo out to radii of 
2.5-5 kpc. The I-front may be preceded by an H 2 photodis¬ 
sociation front driven by LW photons (see e.g. studies by 

). This second radiation 
halo because the rate of 
H 2 formation in the photodissociation region (PDR) exceeds 
the rate at which LW photons are emitted by the star. 

The hard UV spectra of hot high-mass Pop III stars 
increases the thickness of the I-front because of the larger 
range of mean free paths of the ionizing photons in the neu¬ 
tral gas. The outer layers of the front can have tempera¬ 
tures of just a few thousand K and free electron fractions of 


~ 10 %, ideal conditions for the formation of H 2 
phase via the H~ and channels: 

in the gas 

H + e —> H +7 

H" + H —> H 2 + e~ 

( 20 ) 

H + H+ Hj + 7 

Ut + H ->• H 2 + H+ 

( 21 ) 


An H 2 layer may thus form in the outer shell of the front, 
with a molecular mass that greatly exceeds the one in the 
surrounding PDR. In the Hu region itself nearly all H 2 is 
collisionally dissociated by free electrons. Beyond the I-front 
much of the H 2 in the halo that existed before the star was 
born may have been destroyed by the PDR. The total H 2 
column density in the halo may therefore be dominated by 
the H 2 sandwiched between the I-front and the shell of gas 
plowed up by the front after becoming D-type. Since this 
thin layer may govern LW escape from the halo, it is essential 
that it is well resolved in our numerical simulations. 

We show in Figure [2] profiles for the I-front, Hll region 
and PDR for the 60 M@ star in halo 1 at 0.175, 0.35, 0.875 
and 3.5 Myr, the lifetime of the star. At 0.175 Myr it is evi¬ 
dent from the temperature plot that the I-front has become 
D-type: an ionized region with a temperature of 4.5 x 10 4 K 


Ricotti, Gnedin & Shull 2001 2002 


front is often at first confined to the 
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extends out to the position of the front at 1.2 pc, after which 
the temperature drops to 1400 K, rises to 4000 K at 1.5 pc 
and then falls to below 500 K at 1.7 pc. The 4000 K gas is 
the plowed up, shocked material. The temperature falls to 
1400 K in the region between the fully ionized gas and the 
dense shell because of H 2 cooling in the outer layers of the 
I-front. The shock has fully detached from the front and is 
at ~ 1.7 pc. The density spike at 1.3 pc marks the center 
of the plowed up shell, which has ten times the density of 
the ambient halo. H 2 mass fractions reach 10 -4 between the 
I-front and the dense shell and then fall to 10 -7 just beyond 
it, marking the beginning of the PDR. The latter extends 
out to ~ 100 pc, beyond which the H 2 mass fractions grad¬ 
ually rise to ~ 10~ 4 , those that were in the halo prior to the 
birth of the star. Note the extremely low H 2 mass fractions 
within the H 11 region that are due to collisional dissociation. 

At 0.35 Myr, the I-front has advanced to ~ 4 pc, but 
remains D-type and continues to move behind the shock 
front. The structures of the PDR beyond the I-front are 
virtually identical to the ones at 0.175 Myr, as shown in 
the H 2 mass fractions. Note that at each stage the rise of 
the H 2 column density to its peak value coincides with the 
position of the H 2 layer just ahead of the front. The thin 
layer in which most of the H 2 forms has a peak H 2 abundance 
of over 10 —4 , slightly higher than at 0.175 Myr. It is also 
somewhat wider, but the peak density in the layer is roughly 
an order of magnitude lower than at 0.175 Myr. Together, 
these factors lead to the H 2 column density at 0.35 Myr 
being little changed from its value at 0.175 Myr. 

The Hu region remains D-type until about 0.875 Myr, 
when the dense shell is barely visible in the temperature 
profile at ~ 20 pc. The peak H 2 abundance in the shell at 
this point has fallen to 10 -5 because of the smaller densities 
there, resulting in a substantially lower H 2 column density. 
As the halo becomes more and more ionized, more LW pho¬ 
tons pass through the I-front and decrease the H 2 abundance 
in the PDR. This also increases the radius of the PDR, from 
100 pc at 0.375 Myr to 250 pc at 0.875 Myr. The peak den¬ 
sity behind the shock front has fallen to ~ 75 cm -3 , but this 
is still 10 times higher than the density immediately ahead 
of the shock. 

Later on, the I-front becomes R-type, overruns the PDR 
and breaks out of the halo. As the density of the H 2 layer 
falls with the expansion of the I-front, the H 2 column density 
decreases. 

At the end of the lifetime of the star, 3.5 Myr, the whole 
halo is ionized, while the shock front is still moving out¬ 
wards. Nearly all molecular hydrogen is destroyed in the H 11 
region due to collisional dissociations, resulting is a very low 
H 2 column density which is not able to shield against LW 
radiation. 

This case is an example of delayed breakout, in which 
the star eventually ionizes the halo, but only after a large 
fraction of its lifetime. In our grid of models some I-fronts 
loiter at the Stromgren radius until the death of the star. 
For the 9 Mg star in the 1.7 x 10' Mg halo the Stromgren 
radius is below the grid resolution and the Hu region re¬ 
mains hypercompact until the death of the star. In these 
models the halo is essentially undisturbed by UV radiation 
from the star. 


0.0 


0.5 


time (80 M 0 star' 
1.0 1.5 




2.5 


3.0 



Figure 3. Smoothed escape fractions for a 9 Mg star (black 
solid and green dotted lines, overlying) and a 80 Mg star (blue 
solid and red dotted lines) in halo 1 as a function of time. Escape 
fractions with shielding by H 2 only (dotted lines) and by H 2 and 
neutral hydrogen (solid lines) are shown. 


3.2 LW Escape Fractions in the Near-Field Limit 

The escape fraction evolves as the Hu region of the star 
expands in the halo and the H 2 column density and neutral 
H fraction change. In Figure [3] we show how / esc evolves 
over time with H 2 self-shielding alone (dotted lines) and 
with both H 2 and H shielding (solid lines) for the 9 Mg 
and 80 Mg stars in halo 1. With the 9 Mg star, the I- 
front is trapped and ionizing UV photons never break out 
of the halo. With the 80 Mg star, the I-front breaks out 
of the halo and fully ionizes it. Here, the escape fraction is 
essentially zero until the I-front changes to R-type at 0.6 
Myr. Before this transition, the H 2 column density is high 
enough to prevent most LW photons from escaping the halo, 
due to the dense shell driven by the I-front. After the front 
breaks through the dense shell and accelerates down the 
steep density gradient of the halo, the column density of H 2 
rapidly decreases, owing to the rapid decrease in the post¬ 
shock density. Because the density decreases so rapidly with 
radius, most of the H 2 is destroyed by the front well before 
it reaches the virial radius of the halo at a time ~ 0.9 Myr. 

The rapid initial increase in f esc with the 9 Mg star is 
due to the lower luminosity of the star, which creates less H 2 
in the dense shell plowed up by the I-front. The dense shell 
therefore has a lower H 2 column density and hence traps 
LW photons less effectively than in the case of the 80 Mg 
star. However, the smaller UV flux also prevents the I-front 
from overrunning the entire halo, leading to it remaining 
trapped for the lifetime of the star. Consequently, the H 2 
column density in the post-shock shell varies only slightly 
over the lifetime of the star, and the LW escape fraction 
tends towards a roughly constant value. If we account only 
for H 2 self-shielding, then the LW escape fraction at late 
times is around 51%. On the other hand, if we also account 
for absorption in the Lyman series lines of atomic hydrogen, 
we recover instead an escape fraction of around 15% at late 
times. We see therefore that shielding due to neutral H is 
especially important in this case. This is because most of 
the halo is not ionized, and the presence of atomic hydrogen 
reduces / eS c by more than two thirds compared to the value 
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M [Mq] 

9 

15 

25 

40 

60 

80 

120 

halo 1 

45 

13 

34 

58 

72 

78 

83 

halo 2 

34 

8 

< 1 

< 1 

18 

49 

74 

halo 3 

< 1 

10 

34 

3 

< 1 

< 1 

< 1 


Table 3. Near-field escape fractions (in percent) averaged over 
the lifetime of the stars, accounting only for H 2 self-shielding 
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halo 1 
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12 

34 

58 

72 

78 

83 

halo 2 

12 

8 

< 1 

< 1 

18 

49 

74 

halo 3 

< 1 

2 

8 

< 1 

< 1 

< 1 

< 1 


Table 4. Near-field escape fractions (in percent) averaged over 
the lifetime of the stars, accounting for shielding by both H 2 and 
H 
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based on H2 alone. We conclude that shielding by neutral 
hydrogen cannot be neglected for stars with M ~ 10 Mq in 
cosmological halos. 

We list near-field escape fractions that are averaged over 
the lifetime of the star in Tables [3] and [4] Table [3] shows the 
values that we obtain if we only account for H2 self-shielding, 
while Table [4] shows the corresponding values for the case 
where we account for shielding by both H and H2. We see 
that on the whole, we recover similar values in both cases, 
unless the mass of the star is low, or the mass of the halo is 
large. 

Our time-averaged values for the near-field escape frac¬ 
tion are also plotted for all three halos as a function of stellar 
mass in Figure [4] For a given star, f esc falls with increasing 
halo mass because a smaller fraction of the halo becomes 
ionized. In halos 1 and 2, / esc increases with stellar mass, 
because the higher luminosities of the more massive stars 
ionize more of the halo. None of the stars in our study can 
ionize halo 3, which is not only the most massive halo but 
is also the densest (see Figure [l]). All LW photons are also 
trapped, even for the 120 Mq star. With more massive stars 
in halo 3, the Hll region does expand, but the star dies 
before the I-front can reach the virial radius. 

Note that escape fractions for small stars (the 9 Mq star 
in halo, the 9 Mq and 15 Mq stars in halo 2 and the 25 Mq 
star in halo 3) can actually be larger than for more massive 
stars in the same halo. The lower luminosities of these less 
massive stars produce less H2 in the dense shell surrounding 
the Hll region so more LW photons can escape the halo. We 
can compare our numbers with the study by|Kitayama et ah] 
12004). However, we note that they are using an earlier self¬ 
shielding function from DB96 and the overlap in parameter 
space is quite limited. They find the escape fraction for a 
25 Mq star in a 5.5 x 10 5 Mq hal o to be about half the 
value for our 6.9 x 10 5 Mq halo. As Kitayama et al. (20041 
consider their escape fractions to be lower limits, and given 
the differences in halo structure between our models and 
theirs, we conclude that the two approaches yield consistent 
results. 


3.3 Far-Field Limit 

As expected, escape fractions in the far-field limit are gen¬ 
erally much higher than in the near-field. We show time- 
averaged escape fractions with self-shielding due to H2 only 


Figure 4. Time-averaged escape fractions for all our models as 
a function of stellar mass for halo 1 (black thin line), halo 2 (red 
normal line), and halo 3 (blue thick line). Escape fractions with 
just shielding by H 2 (dotted lines) and by H 2 and neutral hy¬ 
drogen (solid lines) are shown, together with some values from 
[Kitayama et ah'](j2004|l (purple stars). 



stellar mass [M s ] 

Figure 5. Time-averaged escape fractions in the far field limit 
for all the models in our study, plotted as a function of stellar 
mass (halo 1: black thin line, halo 2: red normal line, halo 3: blue 
thick line). In the lower panel, escape fractions for self-shielding 
by H 2 alone are shown (dashed lines). In the upper panel, escape 
fractions that account for shielding by H and H 2 (solid lines). 


for all our runs in Table [5] and the upper panel of Figure 
[5] and for shielding by H2 and H in Table [6] and the lower 
panel of Figure [5] respectively. 

When considering only shielding by H 2 , for halos 1 and 
2, the escape fraction is always higher than 90% and nearly 
all LW photons escape the halo. In halo 3, most of the LW 
radiation from the 9 Mq star is trapped in the halo, and 
only a quarter of the photons can escape. The higher mass 
stars all have larger luminosities so their radiation eventually 
escapes halo 3. 

As in the near-field limit, including shielding by atomic 
hydrogen reduces the escape fraction. Less than half of the 
LW photons from the 9 Mq star escape any halo. In par¬ 
ticular the escape fraction in halo 3 is only f eBC « 4%. LW 
escape fractions in halos 1 and 2 for more massive stars are 




























10 Schauer et al. 


M [Mg] 
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98 
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97 
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98 

95 
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97 
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22 

98 

99 

94 

93 

92 
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Table 5. Far-field escape fractions (in percent) averaged over the 
lifetime of the star, accounting only for H 2 self-shielding. 
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15 
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40 

60 

80 
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halo 1 

38 

90 

91 

91 

92 

93 

94 

halo 2 

43 

89 

87 

87 

89 

91 

93 

halo 3 

4 

33 

34 

64 

66 

67 

68 


Table 6. Far-field escape fractions averaged over the lifetime of 
the stars, accounting for shielding by both H 2 and H. 

above 85%, and they rise slightly with halo mass. In halo 3, 
the picture is different: for the 9, 15 and 25 Mg stars, the I- 
front is static and ionizing photons are trapped in the halo. 
The Stromgren radii of the I-fronts of the 15 and 25 Mg 
stars are similar, 0.13 pc and 0.14 pc respectively, which 
explains the step in escape fraction over that mass range. 
With more massive stars the I-front eventually advances be¬ 
yond the Stromgren radius, but the star dies before it can 
reach the virial radius. The final radii of the Hu regions 
increase only slightly with the mass of the star, because al¬ 
though they are more luminous these stars also have shorter 
lives. The radii vary from 12 to 20 pc for 40-120 Mg stars, 
resulting in escape fractions of 64% to 68%. 


4 CONCLUSION 

We have calculated LW escape fractions for 9-120 Mg Pop 
III stars in three different primordial halos ranging in mass 
from 6.9 x 10 s to 1.2 x 10' Mg. We have considered two lim¬ 
iting cases. The near-field limit is relevant for objects close 
to the halo that do not have a significant velocity relative to 
it, so that the LW absorption lines in the rest-frame of the 
object coincide with those in the rest-frame of the halo. In 
this limit, the escape fraction can be approximated by the 
product of the self-shielding function of H2 and the shielding 
function of H2 by neutral hydrogen which have been devel¬ 
oped by WH11. The far-field limit is valid for objects that 
are significantly red-shifted or blue-shifted with respect to 
the halo, usually at cosmological distances. In this limit, we 
can estimate the escape fraction by calculating the ratio of 
the equivalent width of the LW absorption lines to the width 
of the entire LW frequency range. 

In both limits, we consider self-shielding by molecular 
hydrogen alone as well as the combined effects of molecular 
and atomic gas together. We find that it is important to 
consider both types of shielding. 

In the near-field, the escape fraction generally rises with 
increasing mass of the central star. In the 6.9 x 10 s Mg halo, 
it grows from about 10% to about 85% and for the 2.1 x 10 6 
Mg halo to about 75%. None of the stars in our study are 
bright enough to fully ionize the 1.2 x 10' Mg halo and 
escape fractions from this halo are nearly zero in the near¬ 
field. For some stars with masses close to 10 Mg, escape 
fractions can be higher than for slightly more massive stars 
because their lower luminosities cannot build up as large an 


H2 column density, e.g. 12% of the LW photons are able to 
escape the 2.1 x 10 6 Mg halo for a 9 Mg star, but only 8% 
for a 15 Mg star. 

In the far-field, the escape fractions are generally higher. 
Only in the case of the 9 M 0 star do fewer than 50% of its 
LW photons exit all of our halos. In the most massive halo, 
this is also true for the 15 and 25 M 0 stars. For higher stellar 
masses, the escape fraction rises to above 90% for the less 
massive halos and to above 60% for the 1.2 x 10 7 Mg halo. 

If shielding by neutral H is neglected, the escape frac¬ 
tions are overpredicted. Therefore, the additional shielding 
function of neutral hydrogen from WH11 provides an im¬ 
portant improvement to the older model from DB96. This 
is especially severe for stars with masses up to 25 Mg in the 
near-field, where shielding by neutral H reduces the escape 
fraction by as much as a factor of three. In the far-field, less 
than 50% of the LW radiation from a 9 Mg star can exit 
any halo, but if neutral hydrogen shielding is neglected up to 
99% can exit the halo. The 1.2 x 10' M 0 halo is not ionized 
by any of the central stars considered, so its large neutral 
hydrogen column density plays a major role in shielding, re¬ 
ducing the escape fraction from about 90% to about 60% 
for stars with masses of 40 Mg and higher. 

Our models assume ID spherical symmetry. In reality, 
both ionizing and LW photons can break out of the halo 
along lines of sight with lower optical depths, so our escape 
fractions should be taken to be lower limits. We note in 
particular that shielding by neutral H, which has largely 
been ignored in studies until now, can play a large role in 
LW escape fractions, especially those of lower-mass Pop III 
stars. This latter point is important because of the relative 
contributions of lower-mass and higher-mass Pop III stars to 
the LW background. It has recently been found that more 
realistic stellar spectra that contain a lower-mass Pop III 
star component may be less efficient at photodissociating 


H2 in primordial halos than previously thought (Agarwal 
fe Khochfar|2015 1. Shielding by H in the host halos of such 


lower-mass Pop III stars may reinforce this trend by allowing 
fewer LW photons to escape into the cosmos. 

Finally, we note that star formation in the vicinity of a 
halo hosting a low-mass Pop III star may be less influenced 
by the star than by the global LW background because of 
low escape fractions from the halo in the near field. Con¬ 
versely, smaller LW escape fractions from low-mass Pop III 
stars have less of an effect on the LW background on cosmo¬ 
logical scales because escape fractions in the far-field limit 
are generally much larger. 
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